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by : 

L.  A.  Roslund 
and 

N.  L.  Coleburn 


ABSTRACT:  The  stability  to  detonation  and  the  expansion  behavior  of 
reaction  products  were  studied  for  an  aluminized  explosive,  HBX-1. 
The  detonation  velocity-charge  diameter  and  detonation  velocity¬ 
loading  density  relationships  were  determined  by  streak  camera 
techniques  and  electronic  probe  methods .  These  data  gave  reaction 
zone  lengths  which  increase  from  0.07  mm  for  an  initial  explosive 
density  of  1.72  g/cra3  to  1.5  mm  at  1.3  g/cm3 .  With  the  pressure- 
particle  velocity  relation  established  from  338,000  bars  to  15  bars, 
and  conditions  in  the  C-J  state  (Pcj  ■  220. h  kbars  )  specified,  the 
generalized  hydrodynamic  relations  derived  by  Jones  were  used  to 
evaluate  the  constants  in  an  empirical  equation  of  state  and  in 
the  Wilkins  equation  of  state  as  modified  by  Allan  and  Lambourn. 
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I.  INTRODUCTION 


1 .  The  detonation  of  an  explosive  produces  flow  in  the  surroundirg 
medium  which  is  directly  affected  by  the  way  the  available  energy 
is  released  as  the  reaction  products  expand  to  low  pressures.  The 
energy  released  is  from  the  work  integral 

fv 

A  «  I  PdV  ■  Q  -  q  (1) 

Jvo 


where  A  is  the  maximum  available  work  in  an  adiabatic  expansion  of 
the  detonation  products  from  the  initial  explosive  volume  V0,  to  the 
final  specific  volume  v  at  pressure  P.  Q  is  the  heat  of  detonation  and  q 
is  the  heat  retained  by  the  detonation  products  at  V. 

2.  Calculations  of  the  maximum  work  from  Eq.  (1)  are  handicapped 
because  little  is  known  about  the  detonation  reactions  or  the  charac¬ 
teristics  of  the  pressure-volume  expansion  adiabat.  Many  attempts 
have  been  made  to  develop  more  elaborate  equations  of  state  for 

the  detonation  reaction  products.  The  results  of  Deal*  and  con¬ 
siderations  by  Fickett  and  Wood3  seemed  to  support  the  view  that 
for  some  explosives  a  constant  gamma  adiabat  may  adequately  describe 
the  product  gas  expansion  states.  However ,  energy  transfer  computa¬ 
tions  using  the  constant  gamma  assumption  indicate  the  assumption 
is  incorrect,  since  the  results  are  in  considerable  disagreement  with 
experimental  measurements  from  metal  acceleration  experiments. 
Moreover  the  assumption  is  not  expected  to  describe  the  expansion 
behavior  of  burnt  gases  from  explosive  compositions  with  substantial 
after-burning  reactions. 

3.  In  this  work  we  studied  the  hydrodynamic  behavior  of  the  deto¬ 
nation  products  from  an  aluminized  explosive,  HEX-1,  by  making 
measurements  of  the  pressure-particle  velocity  relationship  for  its 
C-J  adiabat.  With  these  data  we  have  described  the  expansion 
behavior  from  33^,000  bars  to  15  bars,  using  a  characteristic  form 
of  the  Wilkins  equation  of  state  (2a)  as  modified  by  Allan  and 
Lambcurn^  and  an  empirical  equation  of  state  (2b) .  These  equations 
are  ; 


•References  are  on  page  13. 
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B  exp  (-kV)  +  WG1V'^1+W)  V  <  Vc  (2a) 

P  = 

g2(v  +  V*)"1  V  >  VC  (2b) 


where  Vq  Is  a  critical  volume  which  specifies  the  properties  of 
the  adiabat.  In  regions  of  high  pressure  with  volumes  less  than 
VC  our  data  are  best  fitted  by  the  first  equation  and  in  regions 
of  expansion  for  volumes  greater  than  Vc  our  data  obey  the  second 
equation.  B,  k,  Gt ,  W,  Gt,  and  V*  are  constants.  The  following 
sections  describe  how  the  experimental  determination  of  the  hydro- 
dynamic  properties  is  used  in  evaluating  the  equation-of- state 
constants.  These  properties  are  the  change  in  the  infinite  diameter 
value  of  the  detonation  velocity  with  loading  density,  the  proper¬ 
ties  of  the  C-J  state,  and  experimental  P-u  data  for  an  expansion 
from  the  C-J  state.  The  expansion  adiabat  derived  from  these  data 
is  compared  with  the  constant  gamma  isentrope  PVY  =  constant,  and 


the  variation  of  y 


with  volume 


is  calculated. 
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II.  EXPERIMENTAL  ARRANGEMENTS  AND  RESULTS 


Charge  Preparation 

4.  The  explosive  configuration  used  in  the  experiments  was  a 
cylindrical  charge  of  pressed  HBX-1,  RDX/TNT/Aluminura/WAX  (40/38.1/ 
17.1/4.8)  by  weight,  which  was  initiated  by  a  plane-wave  expio.ive 
system.  The  pressed  charges*  were  prepared  by  mixing  granular  RDX 

and  TNT  into  a  mixture  of  18  parts  aluminum  and  5  parts  wax  by  weight. 
This  mixture  was  formed  by  stirring  aluminum  and  Stanolind  wax  in 
carbon  tetrachloride  and  then  evaporating  the  solvent.  Seventy 
percent  of  the  HBX-1  by  weight  consisted  of  material  (aluminum,  RDX, 
and  TNT)  with  particle  sizes  between  44  microns  and  180  microns. 
Charges  with  particle  sizes  in  this  range  were  formed  at  the  desired 
charge  density  in  a  hydrostatic  press  and  then  machined  to  the 
appropriate  dimensions . 

Detonation  velocity  Measurements  and  Results 

5.  Detonation  velocities  and  the  critical  diameter  of  HBX-1  were 
measured  as  a  function  of  charge  density  using  smear  samara  and 
"raster  oscilloscope"  pin  techniques.  In  the  smear  camera  measure¬ 
ments,  unconfined  HBX-1  at  loading  densities  ranging  from  about 
1.33  g/cm*  (75%  TMD**)  to  about  1.74  g/cm*  (98%  TMD)  was  detonated 
in  a  "stacked"  cylinder  configuration.  This  configuration  consisted 
of  five  explosive  cylinders  of  different  diameters  but  of  the  same 
loading  density.  The  cylinders  were  stacked  in  order  of  decreasing 
diameters;  5.08-cm,  2.54-cgj,  1.27-cm,  0.635-cm,  and  0.476-cm.  The 
cylinder  heights  corresponding  to  these  diameters  were  10.2-cm, 

10.2-cm,  5.08-cra,  3.81-cm,  and  3.81-cm,  respectively.  The  stacked 
cylinders  were  fired  by  the  detonation  of  a  tetryl  booster  5.08-cm 
diameter  by  2.54-cm  high.  The  tetryl  was  initiated  by  a  5.09-cm 
diameter  pentolite-baratol  plane-wave  booster.  Techniques  described 
in  reference  4  were  used  to  give  narrow,  bright  light  pips  which 
facilitated  the  analysis  of  the  detonation  trace  cn  the  smear 
camera  records.  The  failure  diameter  for  each  charge  density  was 
estimated  by  examining  the  smear  camera  record  for  the  diameter 

at  which  the  detonation  wave  failed  to  propagate  with  a  steady 
velocity.  Single  charges  each  37-cm  long,  5.08-cm  diameter  and 


"*  Cast  charges  also  were  fired  to  compare  their  diameter  effect  with 

pressed  charges  at  the  highest  experimental  density,  about  1.74  g/cm3. 
**  Theoretical  maximum  density,  1.76  g/cm3. 
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2.54-cm  diameter  also  were  fired  simultaneously  with  the  stacked 
cylinder  array  and  the  detonation  velocities  were  measured  using 
electronic  probes  and  a  raster-type  oscilloscope.  In  these  measure¬ 
ments  the  circuit  consisted  of  a  single  strip  of  copper  running  the 
length  of  the  charge,  separated  by  an  air  gap  from  a  series  of 
electrically-charged  precisely  spaced  copper  strips  which  served  as 
probe  contacts  to  be  consecutively  shorted  by  the  passage  of  the 
shock. 

6.  The  camera  data  (which  was  less  precise)  generally  showed 
agreement  with  the  probe  data  to  better  than  0.5%.  The  electronic 
measurements  were  precise  to  better  than  ±0.2%.  With  this  pre¬ 
cision  in  the  velocity  measurements,  the  greatest  source  of  error 
in  the  experiments  probably  is  the  density  variation  within  the 
test  explosive. 

7.  Examination  of  the  detonation  •.-aces  from  the  stacked  cylinder 
array  thowed  that  detonation  did  not  propagate  into  the  0.476-cm 
diameter  charges .  The  critical  minimum  diameter  for  HHX-1  is  about 
0.6  cm.  Prom  the  measurements,  we  determined  a  critical  density 
(between  1.602  g/cm3  and  1.530  g/cms)  below  which  detonation  will 
not  propagate  in  pressed  HBX-1  at  about  0.6-cm  diameter.  We  infer 
that  detonation  cannot  occur  in  smaller  diameters  at  densities 
lower  than  this  critical  density* .  Moreover,  we  found  a  sharply 
different  sensitivity  between  pressed  and  cast  HBX-1.  For  example, 
charges  of  0.63b-cm  diameter  pressed  to  1.72  g/cm3  detonated ; 

oast  charges  of  the  same  density  and  diameter  failed. 

8.  The  detonation  velocities  for  charges  loaded  at  various  densities 
in  the  different  charge  diameters  are  plotted  in  Fig.  1.  Linear 
least  square  lines  are  drawn  through  the  data.  The  effect  of 
charge  diameter  on  the  velocities  is  greatest  at  lower  loading 
densities.  We  extrapolated  the  curves  relating  the  measured  deto¬ 
nation  velocity  to  the  reciprocal  of  the  charge  diameter  for 
different  leading  densities.  For  charge  diameters  >1.27  cm  the 
curves  appear  linear.  The  extrapolation  gave  the  velocities,  !*.,  for 
the  infinite  diameter  or  ideal  curve  which  is  represented  by 

*  -0.063  +  4.305  P  .  (3) 


9.  The  curved  front  theory*  relates  um,  the  reaction  zone  length  a, 
and  the  charge  diameter  d,  to  the  detonation  velocity  D, 

D  -  D„  (1  -  a/d )  (4) 


Assuming  the  equation  only  applies  to  charge  diameters  >1.27  cm, 
we  computed  values  of  the  reaction  zone  lengths  which  are  plotted 
as  functions  of  density  in  Fig.  2.  The  result  is  a  smooth  decrease 
of  the  reaction  zone  length  with  increasing  density.  This  result  is 
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attributed  to  grain  size  reduction  caused  by  the  crushing  action 
of  compaction  to  higher  density.  Smaller  grain  size  means  more 
material  surface  exposed  to  initiation  and  hence  faster  reaction. 

10.  Wood  and  Kirkwood^  give  a  relation  between  reaction  zone  length 
and  curvature  of  the  detonation  front, 

a  =  (R/3.5)  (1  -  D/D.)  ,  (5) 

where  R  is  the  radius  of  curvature  of  the  detonation  front. 

Extensive  studies  of  the  shape  of  the  wavs  by  Cook,  et  al®,  show 
the  front  is  spherical.  Hence,  maximum  curvature,  i.e.,  minimum 
radius,  is  likely  at  the  critical  diameter 

R  *  |  dc  •  (6) 

At  the  critical  diameter,  0.635  cm,  and  critical  density,  about 
1.60  g/cm3,  the  measured  detonation  velocity  is  5234  m/sec,  and 
Do>  ”  6820  m/sec.  These  data  give  a  reaction  zone  length  of  0.21  mm 
as  compared  to  0.40  mm  from  Fig.  2. 


Shock  impedance  Measurements  of  the  Chapman-Jouquet  state 

11.  The  pressure  and  specific  volume  of  the  exDlosion  products  cen¬ 
tered  at  the  C-J  state  were  determined  from  shock  wave  velocity  or  free- 
surface  velocity  measurements,  using  equation-of-state  data  from  several 
sources.*  Shock  waves  were  transmitted  into  several  materials  of  various 
specified  thicknesses  in  contact  with  detonating  HBX-1  and  the  above 
parameters  were  measured  (Table  I) .  To  denote  the  attenuation  and  narrow 
limits  of  the  von  Neumann  spike  region^,  the  data  from  these  experi¬ 
ments  were  linearly  extrapolated  to  values  for  zero  thickness. 

12.  For  these  shock  wave  measurements  the  charge  configuration  consisted 
of  a  12.7-cm  diameter,  6.4-cm  high  cylinder  of  HBX-1,  initiated  on  one 
end  face  by  an  explosive  plane-wave  lens  of  the  same  diameter.  The 
other  face  of  the  cylinder  was  held  in  close  contact  with  the  inert 
material  using  a  thin  silicone-grease  film  to  eliminate  air  gaps. 
Measurements  of  the  free-surface  velocities  and  shock  wave  velocities 
were  made  for  various  plate  and  wedge  thicknesses  of  aluminum,  brass, 
Plexiglas,  and  polyurethane  which  were  shocked  by  detonating  HBX-1  in 
the  above  charge  geometry.  A  reflected-light  smear  camera  technique 
described  previously^  was  used  for  these  measurements.  The  initial 
shock  wave  velocity  in  water  was  measured  from  the  end  face  of  an  HBX-1 
cylinder  5.08-cm  diameter  and  15-cm  long. 

13.  The  C-J  isentrope  behavior  from  about  20  kbars  to  15  bars  was 
determined  from  optical  measurements  of  the  initial  velocities  of 
shock  waves  produced  by  detonating  HBX-1  in  air  and  argon  initially 
at  compressed  or  reduced  pressure  states.  These  experiments  were 
done  in  an  expendable,  gas-tight  chamber  which  could  be  filled 

*See  Table  II  for  equation-of-state  references  for  the  various  materials. 
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with  compressed  gas  at  1000  psi  or  evacuated  to  about  0.01  bar. 

Figure  3  shows  the  experimental  arrangement.  The  chamber  consisted 
of  a  7.62-cm  inside  diameter,  26-cm  long  cast  iron  pipe  tee,  with 
0.5-cm  thick  walls.  The  opposite  ends  of  the  pressure  chamber  tee 
were  fitted  with  2.54-cm  thick  Plexiglas  windows  to  allow  viewing 
of  the  charge  configuration  with  a  rotating-mirror  smear  camera. 

The  bottom  end  of  the  chamber  was  covered  by  a  brass  plate  designed 
to  give  an  air-tight  seal  and  also  to  transfer  detonation  from  a 
detonator  outside  the  chamber  to  a  plane-wave  initiator  inside 
the  chamber . 

14.  The  charges  fired  in  the  chamber  were  5.08-cm  diameter  cylinders, 
16.5-cm  long  and  were  confined  in  copper  cylinders  with  0.64-cm 
thick  walls.  This  confinement  sufficiently  reduces  the  diameter 
effect  on  HBX-1  so  that  the  shock  velocitv  data  are  considered 
applicable  to  an  ideal  charge.  For  viewing  purposes,  one  end 

face  of  the  explosive  protruded  about  1.27  cm  from  the  copper 
cylinder.  The  other  end  face  of  the  HBX-1  was  initiated  by  an 
explosive  train  consisting  of  a  2.54-cm  long,  5.08-cm  diameter 
pentolite  pellet,  a  pentolite-baratol  plane-wave  generator  and  an 
SE-1  detonator. 

15.  The  slit  of  the  smear  camera  was  aligned  to  view  the  shock 
wave  emerging  from  the  protruding  face  of  the  charge  for  a  distance 
of  one  charge  radius,  about  2.54  cm.  Figure  4  shows  the  shock 
wave  trace  emerging  from  a  detonating  charge  in  air  initially  at 
0.00355  g/cm* .  The  camera  writing  speed  was  3.8  mm/ueec. 

The  shock  wave  velocity  was  determined  from  the  distance  vs.  transit 
time  readings  of  the  trace.  Vertical  distance  on  the  film  was 
converted  to  laboratory  coordinates  using  magnification  reference 
markers  positioned  inside  the  test  chamber. 

16.  In  the  region  1  mm  to  b  mm  from  the  charge  surface,  the 
trace  readings  gave  velocities  which  were  characteristically  high 
and  not  reproducible.  This  was  attributed  to  "fog  over"  from 
intense  shock  light  on  the  record,  or  possible  blow-off  of  luminous 
detonating  particles  from  the  charge  surface.  The  region  6  mm  to 
20  mm  from  the  charge  gave  trace  readings  which  were  essentially 
linear.  Slope  measurements  were  made  in  this  region  to  obtain  the 
initial  gas  shock  velocities. 

17.  we  used  the  data  in  refemnce  Li  to  calculate  shock  pressures  jnd 
particle  velocities  cor reSDondina  to  our  measured  air  shock  veloci¬ 
ties*.  Since  our  Masured  shock  velocities  in  argon  are  in  the  range 

of  the  shock  Hugoniot  of  argon  determined  by  Christian,  Duff,  end 
Yarger1*,  we  used  their  relation  for  the  particle  velocities  in 
Table  II.  we  obtained  good  support  for  these  data  from  the  calcula¬ 
tions  by  Bond1*  of  the  shock  Hugoniots  of  argon  for  several  initial 
pressures . 
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HI.  EXPERIMENTAL  DATA  ANALYSIS 


18.  We  obtained  the  data  n  fable  II  from  the  measured  shock  wave 
velocities  and  f  ;ee~surfface  velocities  corresponding  to  zero  thick¬ 
ness  of  the  various  solid  materials  and  their  known  shock  Hugoniots . 
Since  the  pressures  and  particle  velocities  are  continuous  across 
the  reaction  products-material  interface,  the  P-u  data  correspond 

to  points  on  the  detonation  products  Hugoniot  above  the  C-J  state 
and  on  the  expansion  adiabat  below  the  C-J  state.  The  P-u  data 
are  plotted  in  ?ig.  5. 

19.  To  establish  the  C-J  state  we  assumed  that  the  aluminum  and 
wax  in  HBX-1  are  chemically  and  physically  inert  in  the  detonation 
front'4 .  The  detonation  parameters  then  are  functions  of  the  nominal 
density,  pe  of  the  explosive  components  in  the  charge,  i.e.,  the 
ratio  (weight  of  explosive 1/ (volume  available  to  the  explosive). 

For  HBX-1  with  the  experimental  charge  density  of  1.712*  g/cm3,  we 
obtain  =  1.62b  g/cm3. 

20.  .he  C-J  state  was  established  according  to  the  impedance  match 
conditions  by  drawing  the  straight  line  with  slope,  p6D  *  P/u 

th  ough  the  origin  to  intersect  the  curve  given  by  an  exponential 
fit  to  the  P-u  data  above  10  ''.bar.  With  D  *  7307  m/sec  we  obtain 
PCJ  "  220.U  kbare,  u^j  «  18£8  m/sec.  Then  using  the  detonation 
theory,  at  the  C-J  point 

D  -  “CJ  +  C  (7) 


where  C  is  the  velocity  of  sound  at  the  C-J  point .  The  isentropic 
exponent 

YCJ  *  ^P°D~/PCJ^  1 


Also 


Pcj^3 


(9) 


and  for  the  experimental  C-J  point,  *  2.178  g/cm3  and  v  “  2.93b  - 


*This  density  is  typical  for  HBX-1  in  warhead  applications. 
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21.  The  constants  in  the  characteristic  equation  of  state,  Eq.  (2a)  f 
were  evaluated  from  the  experimental  data,  following  the  generalized 
formulation  of  the  detonation  theory  by  Jones18.  Jones  defines  the 
ouantlty  a,  from  the  detonation  products  Shigoniot  by 


3(E  -  Q)p 
~“Tv 


P 

a 


(10) 


where  E  is  the  internal  energy  of  the  detonation  products .  We 
assume  no  radial  expansion  at  the  C-J  plane.  Therefore,  a  is 
evaluated  from  the  infinite  diameter  detonation  velocity-loading 
density  relation,  Eq.  (3)  and  v^j,  the  initial  value  of  the  isentropic 
exponent  in 


..  _  /3  log  P, 

v  "  <rT5gV)s  ‘ 

vs  given  by  Jones,  (reference  15,  Eq.  (14)) 
a  -  (YCJ  +  1)/(1  +  -2 


(id 


(12) 


Pcllowing  Allsn  and  Lambourn's  formulation,  the  constant  W  in  the 
equation  of  state 

I*  ■>  B  exp  (-kV)  +  WGjV'*1  +  W) 


for  V  i  Vc  was  found  from 


W  *  «vCJ/(a  +  1)  (13) 

B  and  Gi  were  determined  from  Bq.  (2a)  and  Eq.  (11)  using  the  C-J 
values  of  P,  V,  and  y.  Differentiation  of  Eq.  (2a)  with  respect 
to  V  gives 

-  H  *  kB  exp  (-kV)  +  W(1  +  WjGjV”*2*"*  .  (14) 

Also,  Eq.  (11)  gives 


dP  _  yP 

‘  3v  T  * 

With  p  *  V”1  and  the  C-J  values  one  then  obtains 

G  »  PCJ(yCJ^CJ  **PcJ  ) 

1  W[ (1  t  W)pCJ  -  kl  » 


(15) 


(16) 
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and 


B  "  Cpcj  -  wg1?cj!1+W’1  exp  <lt/Pcj>  •  (17) 


22.  Then  a  value  of  k  was  chosen  to  give  the  best  fit  to  the 
experimental  P-u  adiabat  in  the  region  V  <  Vg.  To  calculate  the 
P-u  adiabat  from  the  characteristic  equation  we  used  the  Riemann 
relation 


u  -  u 


1 


(18) 


where  C  is  the  velocity  ef  sound;  and  are  the  particle  velocity 
and  density  for  our  known  state,  in  this  case  the  C-J  state.  For 
the  adiabat,  C  is  given  by 


C  * 


(19) 


23.  In  the  region  where  V  >  V~  the  empirical  equation  is 
P  “  Gt/(V+V*).  This  follows  since  in  Fig.  5,  the  experimental 
P-u  data  in  the  lower  region  are  related  best  by 


log  P  *  Mu  +  A  ,  (20) 

where  M  and  A  are  constants  determined  from  the  experimental  data. 
Differentiation  of  Eq.  (18)  and  Eq.  (20)  gives,  with  Eq.  (19), 


»r-s 

(21) 

\dV/  dV 

dP  _  MP 

(22) 

du  dV  „  1 

37  *  Tp”  mp 

(23) 

/l\2/dp\2=  dP 
\M?j  3v 

(24) 

-  -  M2dV 

1 

(25) 

P 


Since  the  reference  state  is  now  the  critical  point  with  the 
characteristic  values  Pc  and  Vc,  we  have 
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(26) 


With 


we  then 


have 


I  .  _L  -  m2<v  -  v  ) 

(27) 

c 

ptv  -  Vc  +  1/(PCM2)]  »  1/M2  . 

(28) 

v*  *  _vc  +  l/(PcM2)  and  Gj  «  1/M2  , 

(29) 

P(V  +  v*)  *  g2 

(30) 

Vc  is  evaluated  by  iteration,  using  Eq.  (18)  to  calculate  the  P-u 
adiabat  from  Equations  (2a)  and  (2b)  and  then  comparing  the 
results  with  the  experimental  data  of  Pig.  5.  Pc  is  calculated 
from  the  characteristic  equation  usina  V  *  Vc.  Equation  (29)  gives 
the  value  of  v*. 
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IV.  DISCUSSION 


24.  The  equation  of  state  constants  for  HBX-1  are  listed  in 
Table  HI  .  The  small  va  .ues  of  a  and  W  are  typical  for  explosive 
compositions  with  large  amounts  of  solid  or  non-ideal  detonation 
products  as  indicated  by  the  relatively  large  values  of  dD/dn0 • 

For  pure  explosives,  a  varies  from  about  0.2  to  0.8,  depending 
on  the  loading  density,  and  dD/dp0  ranges  from  about  3000  to  3500 

JSL/JSL 3 . 

sec  cm 

25.  The  "best  fit"  to  the  experimental  P-u  data  for  V  <  Vc  is 
obtained  with  the  parameter  k  *  6.45  g/cm*;  B  ■  4.187  Mbar  and 

Gt  *  0.00397  Mbar  in  the  characteristic  equation  (Eq.  (2a l  The  experimental 
P-u  relation  for  V  >  Vc  gives  M  -  -10.4589  usec/cm,  G,  ■  9.142  x 
10*  (cm/sec)a.  The  solid  line  in  Fig.  5  is  obtained  when  vc  «  1.15  cnP/g 
then  Pc  ■  3.645  kbar  and  V*  *  1.358  cm3/g. 

26.  If  the  complete  P-u  adiabat  is  calculated  using  the  modified 
Wilkins  equation  for  V  >  Vc,  the  curve  deviates  from  the  experi¬ 
mental  data  as  shown  in  Fig.  5.  Also  plotted  in  Fig.  5  is  the 
constant  y  adiabat  P-u  curve.  The  experimental  P-u  data  from 

20  kbars  to  1  kbar  are  well  below  this  curve.  This  covers  the 
range  of  the  compressed  gas  experiments .  we  note  an  abrupt  upward 
trend  in  the  particle  velocity  data  at  about  10  kbars.  The  trend  is 
considered  real  and  not  a  result  of  experimental  error.  At  points 
below  1  kbar  the  experimental  data  shift  significantly  away  from 
the  constant  y  curve.  These  data  result  from  the  detonation  products 
expanding  into  a  rarefied  atmosphere,  i.a.,  the  initial  conditions 
of  the  reduced  pressure  experiments. 

27.  The  pressure-volume  adiabats  for  the  characteristic  and  empirical 
equation  of  state,  the  modified  Wilkins  equation,  and  the  constant 
gamma  "straight  line"  assumption  are  shown  in  Fig.  6.  The  adiabat 
derived  from  the  experimental  data  dips  below  the  constant  gamma 
curve  near  the  C-J  point  and  then  crosses  over  both  the  constant 
gamma  and  modified  Wilkins  equation  of  state  curves  below  3  kbars. 

The  inflection  in  our  experimental  data  beginning  at  about  3.0  kbars 
produces  the  sharp  break  in  the  adiabat  at  V  ■  1.15  cm3/g.  The 
variation  of  /a  In  p\  with  volume  for  the  modified  Wilkins 

v  "  -\Tnrv)s 

equation  of  state  is  shown  in  Fig.  7.  We  note  that  v  first  in¬ 
creases  from  its  C-J  value  of  2.934,  to  about  3.9,  and  finally 
decreaees  to  1.343.  This  final  value  is  compatible  with  the  value 
of  W  in  the  equation  of  state,  i.e.,  the  limiting  value,  y  ■  1  +  w. 
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It  is  In  the  range  of  expected  values  for  real  gases  at  low 
pressures  as  calculated  from  the  ratio  of  specific  heats ,  Cp/Cy 
about  1.2  to  1.5.  Also  in  Figure  1,  y  derived  from  the  experi¬ 
mental  adiabat  rises  from  its  C-J  value  to  a  maximum  of  5.63 
at  V  *  1.06  cm3/g,  then  decreases  to  5.53  at  V  »*  1.15  cm3/g. 

The  use  of  the  empirical  equation ,  P(V  +  V*)  *  Gj  for 
V  >  1.15  cm3/g  produces  a  discontinuity  in  /d  log  p\  as  in- 

\Hoc  V/S 

dicated  by  the  broken  line  in  Fig.  7. 
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TAbLE  I 

PARTICLE  VELOCITY  VS  THICKNESS  DATA 


Material 

Thickness  (mm) 

u*  ( cm/iisec ) 

Brass 

24.2 

.065 

pc  -  8.40  g/cm3 

24.1 

.066 

1S.0 

.068 

18.0 

.069 

17.0 

.071 

16.1 

.071 

15.0 

.071 

12.7 

.074 

11.0 

.075 

10.5 

.075 

10.0 

.075 

7.62 

.077 

7.11 

.076 

6.22 

.078 

4.83 

.078 

Aluminum 

25.0 

.118 

pQ  *  2.74  g/cm3 

25.0 

.119 

12.3 

.130 

10.0 

.133 

8.0 

.133 

7.0 

.134 

6.4 

.135 

5.0 

.135 

3.2 

.137 

Plexiglas  ... 

12.0 

.212 

p0  ■  1.18  g/cm3 

10.1 

.218 

9.5 

.225 

8.5 

.239 

7.25 

.244 

5.50 

.243 

5.00 

.249 

4.50 

.239 

4.25 

.243 

4.00 

.  243 

3.50 

.239 

3  .00 

.243 

Polyurethane 

14.6 

.237 

»  0.188  g/cm3 

14.4 

.214 

6.45 

.297 

3.38 

.325 

*u  -  Free-iurface  velocity 
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TABLE  II 

PRESSURE -PARTICLE  VELOCITY  MEASUREMENTS 


Material 

Po  (g/cm3 ) 

P0 (bars) 

P (fcbara) 

u (cm/uaec) 

U (cm/u*ec) 

Brass1 

8.40 

3  38 

0.080 

0.503 

Aluminum^? 

2.74 

- - 

270 

0.137 

0.719 

Plexiqlaa18 

1.18 

— 

181.0 

0.243 

0.630 

WateriS 

1.00 

- - 

156 .4 

0.262 

0.597 

Polyurethane^0 

0.188 

— 

37.9 

0.325 

0  .620 

Argon 

7.05 

X 

10-3 

40.8 

16.9 

0.450 

0.534 

N 

4.7 

X 

10-3 

27.0 

11.8 

0.462 

0.54  5 

K 

3.09 

X 

10-3 

17.9 

7.92 

0.466 

0.550 

M 

1.55 

X 

10-3 

9.0 

4.70 

0.510 

0.594 

»• 

1.02 

X 

io-3 

6.0 

3.39 

0.537 

0.621 

Air 

7.10 

X 

io-3 

6.0 

2.54 

0.570 

0.627 

M 

3.55 

X 

io-3 

3.0 

1  .41 

0.602 

0.660 

M 

3.55 

X 

10“3 

3.0 

1.46 

0.613 

0.672 

M 

3.43 

X 

io-3 

3.0 

1.30 

0.587 

0.64  5 

M 

1.27 

X 

10-3 

1.0 

0.67 

0.697 

0.761 

m 

1.27 

X 

10~3 

1.0 

0.67 

0.692 

0.7  56 

m 

8.9 

X 

io-4 

0.7 

0.51 

0.723 

0.794 

m 

6.3 

X 

io-4 

0.5 

0.37 

0.727 

0.793 

M 

6.3 

X 

io-* 

0.5 

0.37 

0.730 

0.796 

m 

3.7 

X 

io-4 

0.3 

0.26 

0.805 

0.875 

M 

3.6 

X 

io-4 

0.3 

0.23 

0.764 

0.832 

*• 

12.6 

X 

10-* 

0.1 

0.101 

0.860 

0.933 

m 

12.6 

X 

10-* 

0.1 

0.101 

0.860 

0.933 

it 

3.86 

X 

10-* 

0.07 

0.069 

0 . 848 

0.920 

m 

8.01 

X 

10~* 

0.07 

0.069 

0.890 

0.965 

n 

1.23 

X 

10~* 

0.01 

0.015 

1 .054 

1.138 

M 

1.23 

X 

10~3 

0.01 

0.015 

1.070 

1.155 

m 

1.23 

X 

10-* 

0.01 

0.016 

1.104 

1.191 

*  Sea  designated  references  for  the  Huyoniot  equation-of-atate 
of  material.  Free-aurface  velocitiea  (»2u)  of  jolid  material* 
were  obtained  by  extrapolation  of  the  data  in  Table  I  to  zero 
material  thickneaa .  Initial  shock  velocitiea  were  meaiiured  for 
water,  air,  and  argon. 
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TABLE  III 

EQUATION-OP- STATE  CONSTANTS  FOR  HBX-1 

1.712  g/c«J 
1.62b  g/cra* 

7.307  mm/u*«c 
2.93b 

220. b  kbar 
l.P?8  nwj/uaac 
7.178  g/cra3 

3p°0 

0.132? 

0.3b32 

1.1?  cm3  /g 

3-6b?  kbar 

6.bc  g/ca* 

b.lfi7  Kbar 

0.00397  Kbar 

9. lb?  x  109  (cm/aac)a 

1.3??  cm3/g 

-10. b''^  uscc/cm 
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FIG.  3  TEST  CHAMBER  FOR  SHOCK  MEASUREMENTS  IN  GASES 
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not 


FIG.  4  SMEAR  CAMERA  RECORD  OF  SHOCK  WAVE  PROPAGATION  FROM  THE  END 
OF  AN  HBX-1  CYLINDER 
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FIG.  5  EXPERIMENTAL  PRESSURE -PARTICLE  VELOCITY  DATA  AND  EQUATION -OF -STATE 
RESULTS  FOR  PRESSED  HBX-1 
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FIG  .  6  PRESSURE  -VOLUME  ADIABATS  FOR  HBX-1 
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'  The  stebility  to  detonation  and  the  expansion  behavior  of  reaction 
products  were  studied  for  an  aluminized  explosive,  HBX-1.  The  deto¬ 
nation  velocity-charge  diameter  and  detonation  velocity-loading 
density  relationships  ware  determined  by  streak  camera  techniques 
end  electronic  probe  methods.  Theee  deta  gave  reaction  zone  lengths 
which  increase  from  0.07  mm  for  an  initial  explosive  density  of 
1.72  g/om3  to  1.5  mm  at  1.3  g/cer3*  with  the  pressure-particle  velocity 
relation  established  from  33?, 000  bars  to  15  bars,  and  conditions  in 
the  C-J  state  (*cj  ■  220.4  kbAr«  )  specified,  the  generalized  hydro- 
dynamic  relations  derived  by  Jones  were  used  to  evaluate  the  constants 
in  an  empirical  equation  of  state  and  in  the  Wilkins  equation  of  state 
as  modified  by  Allan  and  Lambourn. 


“,.1473 


S  -S  3  *  7.  (■  A*"  ? 


UNCLASSIFIED 

Se-  ur»!v 


rac] 


Security  fUsaificstion 


KEY  WORD* 


High  Explosives 
HBX-1 

Aluminized  Explosives 
Reaction  zone  Length 
Failure  Diameter 
Explosive  Effects 
Detonation  Product  Expansion 


DD  ^..1473 


